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We study in detail the model by Isidori and Kamenik that is claimed to explain the top quark 
forward-backward asymmetry at Tevatron, provide GeV-scale dark matter (DM), and possibly im- 
prove the agreement between data and theory in Tevatron W + jj events. We compute the DM 
thermal relic density, the spin-independent DM-nucleon scattering cross section, and the cosmic 
microwave background constraints on both Dirac and Majorana neutralino DM in the parameter 
space that explains the top asymmetry. A stable light neutralino is not allowed unless the local DM 
density is 3-4 times smaller than expected, in which case Dirac DM with mass around 3 GeV may 
be possible, to be tested by the Planck mission. The model predicts a too broad excess in the dijet 
distribution and a strong modification of the missing Et distribution in -I- jj events. 



INTRODUCTION 

There are two questions in modern particle physics 
that have got more attention than any other. The first 
one is the origin of electroweak symmetry breaking. The 
second one is the nature of the cold dark matter (DM) 
of the Universe. The recent CDF measurement of an 
anomalously large tt asymmetry pQ |2] and the obser- 
vation of a peak in the dijet invariant mass distribu- 
tion in association with a W boson [5] have been most 
unexpected and seemingly unrelated to these questions. 
Therefore the proposed explanations to those results of- 
ten involve rather exotic physics meant to explain the 
CDF results alone. 

One very interesting exception in the long list of pro- 
posed models is the one by Isidori and Kamenik [3]. Al- 
though the model was proposed to explain the anomalous 
forward-backward (FB) asymmetry, it involves "usual" 
particles, stops and neutralinos, that may originate from 
supersymmetric models like the next-to-minimal super- 
symmetric standard model (NMSSM). Therefore the 
model provides a DM candidate with mass of order 2- 
3 GeV. (The CoGent experiment recently announced pos- 
sible evidence for light DM [S], but their region of spin- 
independent direct detection cross section is several mag- 
nitudes lower than in the model f^.) Such a light DM 
is very difficult to test even in the most sensitive DM 
direct detection experiments such as XENONIOO and, at 
first sight, seems to be experimentally allowed [HI (see 
[7] for an up-to-date analysis of the XENON results) . In 
addition, the authors claim that the model gives a non- 
resonant contribution to the W + jj channel that may 
improve the agreement between data and theoretical ex- 
pectations. It is intriguing that with a rather minimal- 
istic set of new particles with well-defined couplings and 
masses one can simultaneously explain the tt FB asym- 
metry, the DM of the Universe as well as address the 
observed W + jj excess. 



The aim of this paper is to study the model [4] 
in detail and to work out the related phenomenology. 
We incorporate two versions of the model into the Mi- 
crOMEGAs [SI [3] package allowing the neutralino to be 
either a Dirac or a Majorana fermion. While the col- 
lider phenomenology is not sensitive to the nature of a 
DM fermion, cosmological observables like the DM abun- 
dance and DM annihilation cross section at T = de- 
pend crucially on it. We compute the DM abundance 
and spin-independent direct detection cross section per 
nucleon with MicrOMEGAs by scanning over the par- 
ticle masses and couplings that explain the CDF tt FB 
asymmetry. We find that, indeed, for the model parame- 
ters allowed by the asymmetry, the observed DM thermal 
relic abundance [TU] can be generated and DM as light 
as 2-3 GeV is possible if its couplings to quarks approach 
non-perturbative values. Such a light DM is constrained 
by two observables. First, the observed lack of distortions 
in the Cosmic Microwave Background (CMB) spectrum 
due to DM annihilations at redshifts z ~ 500-1000 puts 
a stringent constraint on GeV-scale DM independently 
of cosmological uncertainties such as the DM halo pro- 
files, densities and distribution [11]. Following Ref. [TT] 
we compute the constraints on {aAv){T = 0) for the an- 
nihilation channel xx ~^ s-^id show that present CMB 
measurements disfavor a light Dirac neutralino as the 
DM. The Majorana neutralino annihilation cross section 
at T = is proportional to the u-quark mass |12j and 
is not constrained. Second, the spin independent DM- 
nucleon scattering cross section is so large that CRESST- 
I 33], TEXONO [HI and XENONIOO [5] experiments 
together exclude the light DM possibility in this model 
entirely. Also, x cannot be a subdominant component of 
DM because this would require increasing its annihilation 
cross section and, consequently, the direct detection cross 
section by the same amount. We conclude that the model 
proposed by Isidori and Kamenik can be a viable model 
for explaining the tt FB asymmetry but the neutralino x 
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cannot be a stable particle and cannot contribute to the 
DM abundance. 

We note that a re-analysis of the experiments [T5] 
presents conservative bounds on detector sensitivities. In 
this case, if, in addition, the local DM density in the 
location of the Earth is a few times smaller than the 
present estimates, a small Dirac neutralino parameter 
space opens up for ~ 3 GeV. The exact measurement 
of the CMB by the Planck mission will be the definitive 
test for this possibility, except when the neutralino is a 
strongly subdominant component (below 30%) of DM. 

In addition, we study whether the model can reproduce 
the dijet excess observed by CDF. We find that the model 
indeed yields an excess of the W + jj signature but the 
invariant mass distribution of the dijet is so broad that it 
cannot explain the CDF observation of a peak, and the 
escaping neutralinos lead to a strong modification of the 
missing Ex distribution. 



THE MODEL 

In addition to standard model (SM) particles, the 
model comprises the color triplet scalar t and neutral 
fermion that are both SU{2) singlets. In the case x° 
is a Dirac fermion, the Lagrangian is 

-CDhac = CsM + {D^t)\DH) - m?ttf + x°(i7M^^)x° 
-mxxV- E (^teix" + H.c.) , 

q—u,c,t 

(1) 

while in the case of Majorana the Lagrangian reads 

q—u,c,t 

(2) 

In numerical computations, the used ranges for param- 
eters were 1.2 < K„ < tt, 4 < Ff < 2tt, 0.1 GeV <m^< 
20 GeV and 194 GeV < < 215 GeV. To avoid bounds 
from flavor physics, we take = 0. The (mf, Yu) pairs 
used were selected from the region compatible with the 
tt-asymmetry at 2cr level in Fig. 1 of [J], extrapolated up 
to the value Yu = tt. We stress that our conclusions do 
not depend on the details of extrapolation. A variation of 
Yt is added for completeness; as long as Yt'^Yu, correct 
ti asymmetry is achieved but the exact value of Yt 
only influences the W + jj and missing Et distribution 
as shown below. Of course, different values of Yt up to 
At: imply different scales for new physics. 



COSMOLOGICAL OBSERVABLES 

The most interesting result of the paper [J is that the 
model simultaneously explains the measured top quark 
FB asymmetry at Tevatron and provides the correct ther- 
mal relic abundance of DM. In order to compute the ob- 
servables related to DM numerically as precisely as possi- 
ble we incorporated both Lagrangians, Eqs. ([T]) and 
into MicrOMEGAs code by using the FeynRules package 
[20] for Mathematica. While for collider observables at 
Tevatron there should not be major differences between 
these two cases, the cosmological observables depend cru- 
cially on the nature of fermionic DM. Because a Majorana 
particle can have only effective four-fermion axial-vector 
interaction |2"lj . both the behavior and the magnitude of 
the DM annihilation cross section will be different in the 
two cases, as well as the DM scattering cross section with 
matter. 

We scan over the model parameters as explained above 
and compute the resulting DM relic abundance and spin- 
independent DM-nucleon scattering cross section. All the 
generated points can explain the CDF tt FB asymmetry. 
Requiring f^DM^^ < 0.1288, which is smaller than the up- 
per limit of the WMAP7 3ct range, we plot in Fig. [1] the 
computed spin-independent nucleon-DM scattering cross 
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FIG. 1: Spin-independent nucleon-DM scattering cross sec- 
tion Gsi (scaled by n^/^DM) vs DM mass m.^ for Dirac and 
Majorana DM. The strength of Yu coupling is represented by 
the color (greyscale) code. XENONIOO i6j, CRESST-I [l3] 
and TEXONO [H] constraints as well as DAMA [TH] best 
fit region are presented. The reanalyzed (r.a.) bounds for 
CRESST-I and TEXONO 115; are shown with dashed hnes. 
The CMB constraints (blue lines) apply only to the Dirac x 
and are given for the case when it is the dominant component 
of DM. 
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section crsi as a function of the DM mass for both 
Dirac and Majorana DM. The magnitude of the coupling 
Yu is represented with the color code. Notice that we do 
allow the model to provide just a subdominant compo- 
nent of the DM of the Universe. As detector sensitivity 
lines are given for the central value of DM abundance, 
the cross section of points is scaled down by the frac- 
tion of X DM in all DM, namely, by ri^/fiDM- How- 
ever, because the annihilation cross section of XX ~^ uu 
and the direct detection cross section of — x^ both 
scale with Yu, the prediction is rather sharp. In the 
same figure we also give the present DM direct detection 
constraints. While XENONIOO excludes DM with mass 
above 5-6 GeV, CRESST-I and TEXONO exclude aU the 
allowed parameter space both for Dirac as well as for Ma- 
jorana DM. Although the sensitivities of CRESST-I and 
TEXONO are not comparable with the XENONIOO one, 
their bounds extend down to 1 GeV because as a target 
CRESST-I used Al and TEXONO used Ge that are much 
lighter than Xe. 

Those model predictions can be readily understood 
with the help of analytic expressions for annihilation 
cross section obtained in [12]. In order to allow the DM 
mass to be in GeV region, the cross section of the annihi- 
lation process xX ~^ and, consequently, the coupling 
Yu, must be large. This explains also why the direct de- 
tection cross sections are large. While the Dirac fermion 
cross section contain both s- and p-wave contributions, 
the Majorana DM annihilation is purely p-wave and, in 
the limit of vanishing u-quark mass, vanishes for vanish- 
ing velocities. Therefore, to obtain a large enough cross 
section for the Majorana neutralino the coupling Y^ must 
be very large. Consequently the spin-independent direct 
detection cross section for the Majorana case is larger 
than in the Dirac case. 

Because the interpretation of results of the direct de- 
tection experiments depends on the local DM density in 
the neighborhood of the Solar System, one can make an 
attempt to save the model by appealing to the conser- 
vative reanalyzed sensitivity bounds for CRESST-I and 
TEXONO given by [TS] and shown in Fig. [l] with dashed 
lines, and reducing the DM density a few times. Indeed, 
cosmological uncertainties in determining the DM halo of 
our galaxy may be larger than expected and the model 
prediction for Dirac neutralino is not too far from the 
CRESST-I bound. However, it was shown in |TT] that 
the present WMAP7 CMB measurements imply strong 
constraints on the annihilation cross section of the GeV- 
scale DM. Those constraints are free of uncertainties of 
DM halo shapes and distributions but apply for the cross 
sections at temperatures T '-^ 0. In our case the CMB 
measurements constrain the Dirac DM case only because 
the Majorana DM annihilation cross section vanishes for 
the reionization temperature T = 10~^ GeV ~ 0. 

We repeated the analyses of Ref. llj for the channel 
XX ~^ The constraints in the {m^,{aAv){T = 0)) 
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FIG. 2: CMB constraints on (ctaw) (scaled by fJ^/f^DM) of 
Dirac DM at T = as a function of DM mass for the 
annihilation channel xX ~^ ww. The coloured (shaded) area is 
the model prediction with the color code for Yu as in Fig. [l] 
Grey lines show the central value of JIdm and the WMAP7 
3(T bounds. 



plane are presented in Fig. [2] together with the area of 
model predictions using the same color code as in Fig. [T] 
The lines in the figure present WMAP7 1 a and 3 a 
bounds and the predicted sensitivity of the Planck mis- 
sion [22 for the central value of JIdm- Because the in- 
tensity of annihilation is proportional to the square of 
density, we scale (ctaw) by the factor of fi^/il^jy;. The 
corresponding CMB constraints are also shown in Fig. [l] 
for X as the dominant component of DM. Consequently, 
if the detector bounds could be reduced at least by 3-4 
times by assuming that the local DM density is smaller 
than expected (in effect, detection sensitivity lines shift 
up by the same factor), at 2 cr level a very small pa- 
rameter space at = 3.3-3.6 GeV for Dirac DM opens 
up. However, the variability of the local density of DM 
is unknown experimentally and very poorly estimated by 
computer simulations at the spatial scales of the Solar 
system [17l [18] . Thus we are unable to give a quantita- 
tive estimation for the probability of having local density 
of DM that is a few times smaller, but such a reduction 
is highly unlikely (e.g. [19)). This parameter region will 
be definitively tested by the Planck mission, except when 
neutralino DM constitutes less than 30% of all DM. 

To allow Majorana DM to be viable, the local DM 
density must be reduced by 2 orders of magnitude, which 
is not realistic. 



W + DUET DISTRIBUTION AT TEVATRON 

The CDF collaboration recently published data on the 
invariant mass distribution of jet pairs produced in as- 
sociation with a W, finding an excess of events around 
150 GeV after background subtraction [3]. The DO col- 
laboration, on the other hand, has seen no such signal 
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As noted in [4] , the model contributes to this signature 
via stop pair production, 



xu) w+b)) . 



(3) 



We simulate this channel as well as the Standard Model 
diboson contribution for two benchmark points with Yt — 
3 and 4, respectively, fixing in both cases y„ = 1.5 and 
mj = 205 GeV to reproduce the tt FB asymmetry and 
— 2 GeV (the results are weakly dependent on the 
precise value of m^). We use MadGraph 5 [24j for parton- 
level event generation, Pythia 6 [ 25] for hadronization 
and PGS 4 width the CDF parameter set for an esti- 
mate of detector effects. We apply the kinematical cuts 
described in T at the detector level. 

The distribution of the dijet mass rrijj for the Standard 
Model and the two benchmark scenarios are shown in 
Fig. |3] The model indeed leads to an excess in the same 
region as seen in the CDF analysis. However, while we 
do not directly compare the simulation to the data, in 
view of our imperfect detector simulation, lack of next- 
to-leading order and jet energy scale corrections, it is 
obvious that the distribution is much broader than the 
CDF bump and that the excess thus extends to much 
lower and higher values of rrijj . 

In addition, the two neutralinos present in the final 
state lead to a significant increase in the missing energy 
per event. We show the missing Et distribution, after ap- 
plying the same cuts as for Fig. [3j for the diboson contri- 
bution and the two benchmark points in Fig.|4] While the 
corresponding background-subtracted plot from the ex- 
periment is not public to our knowledge, it should clearly 
show this huge excess, to be tested at the LHC as well 
(the CMS detector resolution for missing Et is around 
10 GeV [IS]). 

We conclude that, while the model might slightly re- 
duce the tension between theory and data for the nijj 
distribution, the strong modification of the missing Eq^ 
distribution excludes this possibility. 
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FIG. 3: Dijet invariant mass rrijj distribution in the two 
benchmark scenarios (blue, light grey) and the Standard 
Model diboson contribution (red, dark grey). 
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FIG. 4: Missing Et distribution in the two benchmark sce- 
narios (blue, light grey) and the Standard Model diboson con- 
tribution (red, dark grey). 



CONCLUSIONS 

We have performed a detailed study of the Isidori and 
Kamenik model [Ij. We find that the DM direct detec- 
tion experiments and the CMB measurement by WMAP7 
exclude the possibility that the Dirac or Majorana neu- 
tralino could contribute to the DM relic abundance. If 
the local DM density in the Solar System is reduced 3-4 
times (a rather unlikely possibility), a small parameter 
region for — 3.3-3.6 GeV may open for DM at 2 ct 
level. This possibility will be tested by the Planck mis- 
sion, except when neutralino DM makes up less than 30% 
of all DM. Alternatively the model must be modified so 
that the neutralino x decays before or after nucleosyn- 



thesis. In this case the model still provides a viable sce- 
nario to explain the ti FB asymmetry but not DM. The 
model cannot explain the CDF dijet anomaly due to a 
too broad dijet invariant mass distribution and significant 
additional missing transverse energy in -I- jj events. 
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